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This  report  illustrates  that  cogeneration  plant  operation  can  be  optimized  through  the  use 
of  a  hand-held  programmable  calculator.  The  cogeneration  system  consisting  of  four  high-pressure 
steam-generating  units  and  three  turbine  generators  at  U.  S.  Naval  Submarine  Base,  New  London, 
CT  is  analyzed.  A  computer  program  written  for  use  on  a  TI-59  programmable  calculator  is  there¬ 
fore  developed  to  provide  the  personnel  of  the  power  plant  a  rapid  determination  of  the  optimum 
operating  settings  for  the  turbine/generator  sets  at  any  given  time.  Economic  figures  are  presented 
to  emphasize  that  the  implementation  of  this  program  would  result  in  significant  improvements  in 
economics  and  efficiency  of  the  cogeneration  process  as  the  fuel  price  is  expected  to  increase  more 
rapidly  than  the  utility  rate.  Although  this  computer  program  is  initially  developed  for  the  appli¬ 
cation  at  SUBASE,  New  London,  CT,  its  applicability  to  other  similar  cogeneration  systems  can 
be  achieved  by  simple  modification  of  the  site-specific,  sub-routine  functions. 
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INTRODUCTION 


The  importance  of  cogeneration  has  been  greatly  emphasized  as  the 
cost  of  energy  has  increased.  The  cost  of  purchased  electrical  power 
has,  of  course,  risen  with  fuel  prices.  The  increases  in  investment 
cost  for  power-producing  equipment  have  forced  utility  companies  to 
develop  new  rate  structures  which  emphasize  the  reduction  of  reserve 
capacity  by  severely  penalizing  peak  usage.  All  this  has  resulted  in 
careful  reviews  of  internal  power-generating  capabilities. 

Currently,  the  U.S.  Navy  has  approximately  12  operational  steam- 
electric  power  plants  with  extraction  turbines.  They  are  operating  in 
parallel  with  the  utility  company  serving  the  area.  These  plants  are 
representative  of  numerous  small  industrial  turbine  facilities  that  are 
capable  of  practicing  cogeneration.  A  major  concern  in  the  cogenerating 
power  plant  has  been  that  the  operation  of  the  extraction  steam  turbine 
electric  plants,  generally  based  on  operator's  judgments,  may  fail  to 
achieve  maximum  economy  in  power  usage  at  Naval  installations.  At  some 
arbitrary  control  settings,  or  in  response  to  a  demand  level  at  which 
peak  shaving  is  initiated,  some  electrical  power  for  the  activity  is 
"bought"  from  a  utility  company  and  some  power  is  "made"  onsite.  This 
combination  of  purchased  versus  self-generated  power  may  not  result  in 
the  least  cost  power  mix  as  power  demand  and  steam  loads  vary. 

The  need  for  a  means  of  determining  the  optimum  economic  mix  of 
self-generated  and  purchased  power  in  a  plant  that  cogenerates  steam  and 
electricity  has  long  been  recognized  by  the  U.S.  Navy.  In  view  of  the 
possible  near-term  benefits  which  can  be  derived,  the  Civil  Engineering 
Laboratory  (CEL)  was  funded  to  develop  a  computer  program.  This  computer 
program  will  provide  the  plant  operator  with  a  method  of  rapid  determina¬ 
tion  of  the  optimum  make/buy  operating  decisions  for  the  Navy's  existing 
automatic  extraction  steam  turbine/generator  cogenerating  plants.  The 
objectives  of  the  study  are  to: 


1.  Select  a  potential  Naval  demonstration  site  (U.S.  Naval  Submarine 
Base  (SUBASE)  at  New  London,  Conn,  was  selected). 

2.  Develop  an  algorithm  for  determining  the  cost  of  self-generated 
electricity. 

3.  Formulate  the  site-specific  functions  for  determining  the  cost 
of  purchased  power. 

4.  Develop  a  computer  program,  based  upon  the  above  results,  for 
use  on  the  Texas  Instruments  TI-59  programmable  calculator  with  printer 
for  determining  the  most  economical  mix  of  self-generated  and  purchased 
electricity . 

5.  Provide  user  with  instruction  on  the  use  of  the  computer  program. 

Although  the  computer  program  presented  in  the  report  is  developed 
for  the  application  at  SUBASE,  New  London,  Conn.,  it  can  be  modified 
according  to  site-specific  information,  such  as  the  system  characteristics 
and  utility  rate  structure,  for  the  application  at  other  Navy  activities 
with  installed  steam  turbine/generator  equipment. 

The  assistance  and  cooperation  of  personnel  at  the  Public  Works 
Office  at  SUBASE,  New  London,  Conn,  is  greatly  appreciated,  particularly 
that  of  CDR  James  C.  Hay,  LCDR  Alan  H.  Burkett,  and  Mr.  Max  C.  Browning. 
Special  recognition  is  given  to  LCDR  Alan  Burkett,  CEC,  USN,  for  his 
assistance  in  optimizing  the  computer  programs  for  both  user  convenience 
and  run  speed. 


SYSTEM  DESCRIPTION 

The  cogeneration  system  at  SUBASE,  New  London,  Conn,  is  depicted  in 
Figure  1.  It  consists  of  four  high-pressure  steam-generating  boilers 
and  three  turbine  generators.  All  of  the  boilers,  each  rated  at 
76,500  pounds  of  steam  per  hour,  are  fueled  with  No.  6  oil  and  supply 
steam  to  600-psig  headers  which  are  interconnected  by  crossties.  Steam 
in  the  600-psig  header  is  usually  heated  to  about  700°F.  Boiler  no.  IS 


is  a  newer  boiler  than  the  other  three.  The  overall  efficiency  of  these 
boilers  is  assumed  to  be  72%.  The  makeup  for  all  heat  balances  is 
assumed  to  lie  between  17  and  19%  of  steam  generated  by  the  boilers,  and 
3%  of  this  is  assumed  to  be  blowdown.  The  majority  of  the  losses  repre¬ 
sented  by  the  difference  between  the  makeup  and  the  blowdown  is  assumed 
to  be  lost  in  the  export  system. 

All  three  turbines  are  condensing  units  with  automatic  extraction. 
Manufacture  information  for  the  turbine/generator  sets  is  listed  in 
Table  1.  Turbines  no.  3  (TG3)  and  no.  5  (TG5)  receive  throttle  steam  of 
600  psi,  700°F  from  the  boilers  and  both  extract  200-psig,  500°F  steam 
for  the  header  to  turbine  no.  4  (TG4),  to  the  steam-driven  power  plant 
auxiliaries,  and  to  the  200-psig  export  supply.  The  5-psig  steam  is 
supplied  by  the  bleed  from  the  5-psig  extraction  ports  of  turbines  no.  3 
and  no.  4.  A  pressure-reducing  valve  (PRV)  and  a  desuperheater  (DSPH) 
station  are  piped  between  the  600-  and  200-psig  headers  for  the  direct- 
pressure  reduction  purpose.  There  is  normally  sufficient  extraction 
steam  capacity  in  the  turbines  so  that  the  PRV  can  remain  closed.  If, 
for  any  reason,  the  steam  pressure  in  the  200-psig  header  drops  below 
the  limit,  the  PRV  will  open  to  maintain  the  header  pressure.  Performance 
curves  for  the  three  turbines  are  shown  in  Figures  2  to  4  from  which  the 
overall  turbogenerator  efficiencies  are  estimated  and  listed  in  Table  1. 

The  self-generated  power  onsite  and  the  purchased  power  from  the 
Groton  City  Utility  Company  are  connected  in  parallel  to  supply  the 
switchgear  in  the  power  plant. 


SELF-GENERATION  ANALYSIS 

The  algorithms  to  be  employed  to  calculate  the  cost  of  electrical 
power  generated  by  the  cogenerators  onsite  are  presented  in  this  section. 

Considering  each  turbine/generator  set  shown  in  Figure  1  as  a 
control  volume,  equations  which  represent  the  conservations  of  mass  and 
energy  principles  can  be  written  as  follows: 
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Continuity  equation: 
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First  law  of  thermodynamics: 
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The  explanation  of  the  symbols  used  in  the  above-mentioned  equations 
and  in  Figure  1  is  presented  in  the  Nomenclature  List. 

Combining  these  equations  gives: 


+ 


^M3H  +  M3L  +  M3E^  H3T  +  (M4L  +  M4E^  H4T 

+  (M5H  +  M5E^  H5T  "  (M3H  x  H3H) 

(M3L  x  h3L)  -  (M3E  x  H3£) 

'  <M4L  X  H4L>  '  (M4E  X  H4E) 

'  <M5H  x  »5H>  -  <M5E  X  (7) 
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Assuming  that  the  piping  losses  jic  negligible,  then 


H3T 

H5T 

Hs  . 

(8) 

H3H  = 

H5H  = 

H4T  =  HH 

(9) 

Hu  = 

hm  - 

H1 

(10) 

Subst i tut  ion 

ot  t  he 

se  irlji  i.iii'  .ilong  with  the  following  identities 
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From  the  heat  balance  of  this  system,  it  has  been  found*  that  the 

quantity  of  the  5-psig  steam,  to  be  supplied  to  the  deaerator,  is  about 

the  same  as  that  to  be  exhausted  from  the  three  steam-driven  auxiliaries 

as  shown  in  Figure  1.  Therefore,  it  is  reasonable  to  assume  that  + 

is  equal  to  M^,  the  low-pressure  export  steam. 

Assume  that  the  amount  of  steam  consumed  for  steam-driven  auxiliaries 

between  the  200-psig  and  the  5-psig  headers  is  negligible  compared  with 

high-pressure  export  steam  The  term  (M^  +  ^5^  -  M^)  Equation 

11  will  be  equal  to  the  high-pressure  export  steam  supply  Mu.  In  addition, 

n 

it  is  reasonable  to  assume  that 


*Carlson  &  Sweatt  Engineers.  "Report  on  economics  of  electrical  power 
generation,"  New  York,  N.Y. ,  Aug  1972.  (Prepared  for  Naval  Submarine 
Base,  New  London,  Conn.) 
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HS  '  H3E 


* 


4E 


HS  '  H5E 


With  the  above-mentioned  assumptions,  Equation  11  can  be  rewritten 
as 

hh  <hs  -  hh>  *  "l  (hs  -  V 

*  «3E  ♦  M«  +  "sE*  <"s  -  "e>  (12) 

Where  H„  is  defined  as  the  enthalpy  of  the  mixed  exhaust  steam,  and 
h 

M0„  +  M/t,  +  Mc  is  the  total  condensing  steam  flow  of  this  system  with 
3h  4h  jL 

the  limit  value  of  M_..  +  M...  +  Mru. 

3M  4M  5M 

Up  to  this  point,  Equation  12  has  been  derived  for  the  total  thermal 
energy  to  be  converted  into  the  predetermined  electrical  power, 

+  W  while  meeting  the  steam  load  demand  at  any  given  time. 

In  order  to  compute  the  cost  for  the  electricity  generated  onsite, 
the  turbine  heat  rate  (T  )  is  determined  and  used  as  follows: 

Total  energy  added  to  turbines: 

(mjt  +  m4t)  hs 

The  thermal  energy  which  is  recoverable  from  the  turbine  system: 

MH  X  HH  +  ML  X  HL  +  ^M3E  +  M4E  +  M5E^  x  HC 

Where  H^,  is  the  average  enthalpy  value  of  the  condensate  from 
turbines . 

The  turbine  heat  rate  (T™)  is  the  net  heat  input  to  the  turbines 
per  unit  of  electrical  generated,  which  can  be  written  as: 

THR  =  ^M3T  +  M5T^  HS  "  ^MH  X  HH^  "  ^ML  X 

'  (M3E  +  M4E  +  M5E)  HC]/(W3  +  W4  +  V  (13) 
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Substitution  of  Equation  (12)  into  Equation  (13)  leads  to 


UK 


I  <w  ,/ 1-: )  +  ( \< .  /!■:,  ,.)  +  <w  / 1:  , )  |  i 

5)1  -t  -4  1  )  I  1 


+  (M  jl;  +  Mti.  +  M5[;)  (II,-  -  H(J 


+  U’  +  W  ,  I 

-4  J 


(14) 


Where  the  turbine  heat  rate,  T™ ,  in  unit  of  Btu  per  kW-hr,  is  the 

ni\ 

instantaneous  heat  rate  derived  from  a  particular  heat  balance  of  the 

overall  cogenerating  system  at  a  fixed-load  situation. 

Examination  of  Equation  14  indicates  that  it  will  be  more  efficient 

to  generate  electricity  onsite  when  the  system  is  operating  at  the 

minimum  allowable  condensing  level,  i.e.,  when  (Mor  +  M,„  +  Mc„)  is 

Jh  4h  Dh 

minimum.  It  is,  therefore,  the  objective  of  this  study  to  seek  an 
optimal  balance  between  the  power  demand  set  point  and  the  operation  of 
the  three  condensing  turbines. 

With  the  information  available  for  the  heating  value  of  fuel,  fuel 
price  and  boiler  efficiency,  the  cost  of  self-generation  of  W^  (=  W^  + 

W4  +  W^)  electrical  power  at  any  expected  export  steam  loads,  and 
can  be  expressed  as  follows: 

Unit  cost  of  self-generation,  C/kW-hr  is 


JUS 


=  thk  x  V(Hv  *  eb> 


(15) 


Total  cost  of  self-generated  electricity,  W^  +  W^  +  W^ ,  over  the  time 
period  of  T  hours  (which  may  be  a  period  of  one  shift  or  less)  is 


JTS 


=  CUS  X  (W3  +  W4  +  W5)  x  T 


(16) 


PURCHASED  POWER  ANALYSIS 


The  SUBASE  at  New  London,  Conn,  purchases  part  of  its  electrical 
power  from  the  Department  of  Utilities  of  the  City  of  Groton.  The  rate 
structure,  presented  in  Table  2,  includes  two  charges  during  the  30-day 
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billing  period:  the  demand  charge  and  the  energy  charge.  The  energy 
charge  is  simply  a  predetermined  price  per  kW-hr  used.  The  price  can 
vary  with  the  amount  purchased,  but  basically  it  is  a  straightforward 
cost.  The  demand  charge  is  currently  established  as  a  predetermined 
price  times  the  highest  peak  value  created  over  the  last  11  months.  The 
demand  charge  represents  a  very  significant  part  (about  25%)  of  the 
total  utility  bill.  In  such  a  structure,  the  demand  charge  becomes 
almost  like  a  fixed  cost,  which  has  to  be  paid  whether  it  is  utiliaed 
or  not.  Each  bill  for  the  current  month  will  be  increased  or  decreased 
by  an  amount  equal  to  the  fuel  adjustment  factor  (F.)  times  the  total 

Pi 

energy  (kW-hr)  used  in  the  current  billing  month. 

The  cost  of  the  purchased  electricity  over  the  time  period,  T, 
during  the  billing  month  is  calculated  according  to  the  rate  structure 
listed  in  Table  2  as  follows 

(i)  The  demand  charge,  $/kW-hr: 

4.72/(30  x  24) 

(ii)  The  energy  charge,  $: 

(0.0343  +  Fa)  +  (0.0269  +  FA)  R2  +  (0.0211  +  FA)  R3 

The  values  of  coefficients  R^,  R2  and  R^  are  dependent  upon  the 
values  of  the  cumulative  energy  consumption  at  the  beginning  (K^),  and 
that  at  the  end  (K2)  of  the  time  period,  T.  The  relation  between  and 
can  be  presented  as 

K2  -  Kj  =  Kp  x  T 

Where  Kp  represents  the  electrical  power  to  be  purchased  during 
time  period,  T,  from  the  utility  company. 

These  coefficients  Rj ,  R2  and  are  to  be  determined  according  to 
the  following  conditions: 
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a 


A 


(1)  K2  <  50,000 

Rj  =  Kp  X  T 

r2  =  s  =  0 

(2)  50,000  <  K2  <  300  x  Du 

(a)  Kj  <  50,000 

Rj  =  50,000  -  Kj 

R2  =  K2  -  50,000 

r3  =  0 

(b)  Kj  >  50,000 

R2  =  Kp  x  T 


(3)  300  x  Du  <  K2 

(a)  Kj  <  50,000 

Rx  =  50,000  -  Kx 

R2  =  300  x  Dy  -  50,000 

R3  =  K2  -  300  x  Dy 

(b)  50,000  <  K1  <  300  x  Du 
R1  =  0 

R2  =  300  x  Du  -  Kj 
R3  =  K2  -  300  x  Dy 

(c)  300  x  Dy  <  Kj 

R1  =  R2  =  0 
R3  =  Kp  x  T 


Thus,  the  total  cost  (C^,p  in  dollars)  for  the  purchased  utility 
(Kp)  over  the  period,  T,  is 
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(17) 


CTp  =  (0.0343  +  F  )  Rj  +  (0.0269  +  Fft)  H2 

+  (0.0211  +  Fa)  R3  +  (4.72/(30  x  24)]  x  x  T 

And  the  unit  cost  of  utility  (C^p  in  C/kW-hr)  for  the  purchased  power 
(Kp)  is 

CUP  =  Ictp/(kp  x  t)1  x  100  (18) 


OPTIMUM  MAKE/BUY  DECISION  ANALYSIS 


When  there  is  a  combination  of  self-generation  and  purchase  of 
electrical  power,  the  total  unit  cost  (in  cents  per  kW-hr)  can  be 
expressed  as: 


CTS  *  CTP 
U  D  x  T 


(19) 


Where  C^, and  C^,p,  the  total  costs  of  self-generation  and  purchase 
power,  respectively,  have  been  derived  in  the  preceding  sections.  It  is 
evident  from  Equations  16  and  17  that  at  the  specified  demands  of  power 
and  steam,  both  C^g  and  C^,p  are  a  function  of  the  total  self-generation 
onsite,  W^  (=  W3  +  W^  +  W^ ) .  Thus,  Equation  19  can  be  written  mathe¬ 
matically  as: 

Cy  (WT)  =  CTS(WT)  +  CTp  (WT)  (19-1) 


A  realistic  problem  of  this  combination,  Make  versus  Buy,  which  has 
been  encountered  by  plant  operators,  is  how  much  power  ought  to  be 
generated  onsite  in  order  to  meet  the  requirements  at  the  lowest  cost. 

In  order  to  solve  this  optimization  problem,  which  is  the  objective  of 
this  study,  it  is  necessary  to  select  the  optimum  value  of  self- 
generation,  W,p,  to  minimize  the  total  unit  cost  derived  from  Equation  19. 
A  necessary  condition  for  a  minimum  value  of  to  exist  is  that: 
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=  0 
T 


(19-2) 


This  is  done  here  numerically,  using  a  hand-held  programmable 
calculator,  in  accordance  with  the  following  operational  limits  for  the 
system : 

•  Maximum  extraction  flow  for  each  steam  turbine 

•  Minimum  condensing  flow  for  each  steam  turbine 

•  Limits  for  turbine/generator  power  output 

It  should  be  noted  that  the  solution  of  determined  from  Equation  19 
may  not  be  minimum  where  Equation  19-2  holds,  but  an  optimum  one  with 
which  the  above-mentioned  constraints  comply. 

The  numerical  procedure  for  determining  the  optimum  value  of  unit 
cost,  Cjj ,  starts  with  the  assumption  that  maximum  possible  electrical 
power  will  be  generated  onsite  and  any  additional  electricity  will  be 
purchased  in  order  to  meet  the  power  demand  requirement  over  the  time 
period,  T.  The  resultant  unit  cost,  C^,  is  then  compared  with  the  new 
value  resulting  from  the  assumption  of  reducing  the  self-generation 
capacity  by  a  predetermined  decrement  and  the  increasing  purchased 
electrical  power.  This  iteration  is  continued,  i.e.,  the  assumption  of 
generating  less  electricity  onsite  and  purchasing  more  from  the  utility, 
until  the  optimum  (or  minimum)  value  of  unit  cost  is  found. 


COMPUTER  PROGRAMS 

Program  Description 

The  computer  programs  are  developed  according  to  the  analysis 
presented  in  the  previous  section  and  are  designed  for  use  on  a  Texas 
Instruments  TI-59  programmable  calculator  with  a  PC-100C  printer.  The 
program  can  be  utilized  by  the  operators  of  the  power  plant  at  any  time 
to  determine  the  optimum  mix  of  self-generated  and  purchased  electrical 
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power  at  SUBASE,  New  London,  Conn.  The  results  of  this  program  will 
provide  a  guideline  which  will  enable  the  operators  to  operate  the  power 
plant  in  the  most  economical  way. 

A  flow  chart  describing  the  logic  of  the  computer  program  is  shown 
in  Figure  5.  The  listing  of  the  computer  program  suitable  for  use  on 
the  Texas  Instruments  TI-59  programmable  printing  calculator  appears  in 
Table  3.  The  first  two  columns  of  the  program  listed  in  Table  3  are  for 
the  program  location  and  key  code,  respectively.  The  group  of  three 
digits  shows  the  location  in  program  memory  of  each  instruction.  This 
not  only  allows  the  user  to  keep  track  of  instructions,  but  also  tells 
the  calculator  the  order  in  which  to  complete  the  instructions.  Since 
the  calculator  can  only  understand  numbers,  each  key  symbol  on  the 
keyboard  is  assigned  a  two-digit  code  number  known  as  a  key  code  corre¬ 
sponding  to  the  instruction  stored  in  the  program  memory  location. 

There  are  two  groups  of  information  required  for  this  program  to 
run.  The  first  group  of  information  is  the  characteristics  of  the 
onsite  generating  system  which  can  be  named  "System  Information"  including 
the  efficiencies  of  turbine/generator  sets  and  boilers,  the  thermodynamic 
properties  at  the  points  of  interest,  and  the  constraints  of  the  turbine/ 
generator  sets.  All  the  systems'  information  is  treated  as  constants 
and  is  stored  at  data  registers  0  to  17.  These  values  are  listed  in 
Table  4.  If  any  change  is  necessary  due  to  equipment  modifications,  it 
is  relatively  easy  to  revise  the  required  values  with  a  basic  under¬ 
standing  of  programming  on  the  TI-59  calculator.  It  is  worthwhile  to 
mention  that  the  constraints  of  the  turbine/generator  sets  are  stored 
with  the  assumption  that  all  three  sets  are  on-line  for  operation. 

However,  the  user  still  has  the  options  to  make  his  choice  of  any  other 
combination  of  turbine/generator  sets  in  the  process  of  determining  the 
optimum  mix  of  self-generated  and  purchased  electrical  power.  To  assume 
putting  an  out-of-service  turbine/generator  set  back  in  service,  the 
user  has  to  take  some  simple  steps  to  restore  into  the  data  registers 
the  constraints  of  the  equipment  before  the  user  starts  to  run  the 
program.  Details  of  this  procedure  will  be  described  later. 
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»i  Information," 


The  second  group  of  information,  called  "Operat 
consists  of  ten  variables.  These  all  have  to  be  keyed  into  the  calculator 
by  the  user  when  the  program  is  run  for  the  first  time.  For  subsequent 
runs,  only  those  variables  which  are  different  from  the  previous  ones 
need  to  be  keyed  in  again  with  the  exception  that  the  values  of  W^, 
and  W,_  have  to  be  reentered  for  every  run.  The  ten  variables  and  their 
data  registers  are  listed  as  follows: 


*  W3 

*  W4 

*  W5 

•Mh 

*  ml 

*  D 

*  K 


1 


•PB 

•  T 


R20,  highest  selected  electric  output  (kW)  of  TG3 

R21,  highest  selected  electric  output  (kW)  of  TG4 

R22,  highest  selected  electric  output  (kW)  of  TG5 

R2^,  high  pressure  (200  psig)  export  steam  (lb/hr) 

R2^,  low  pressure  (5  psig)  export  steam  (lb/hr) 

R2^ ,  electric  demand  (kW)  during  the  time  period,  T 
R^,  cumulative  purchased  energy  (kW-hr)  at  the  beginning 
of  each  run 

R2^,  fuel  adjustment  factor  ($/kW-hr) 

R2g,  fuel  price  (C/gal.) 

R2g,  time  period  of  interest  (hr)  for  each  run 


All  these  ten  inputs  are  stored  automatically  by  the  program  into 
data  registers  20  to  29  (^20"P29^'  ^rst  t^ree  inputs  are  the 

selected  generator  outputs  of  turbine/generator  sets,  which  are  normally 
the  rated  capacities,  with  the  condition  that  the  total  inputs  of  W^, 
and  W,.  at  any  time  must  be  either  equal  to  or  less  than  the  input  of 
electrical  demand  D. 


Operating  Instructions 

A  detailed  explanation  of  how  to  run  this  program  is  presented  in 
this  section.  A  step-by-step  user  instruction  to  run  the  program  is 
shown  in  Table  5.  The  explanation  of  each  step  is  given  below. 

j 

■i 

) 

:1 

i 
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Slide  the  ON/OFF  switches  on  the  T 1—39  calculator  and  the 
printer  to  the  ON  positions.  A  single  zero  should  be  seen  in 
the  calculator  display. 


2  Read  magnetic  cards.  The  two  magnetic  cards  used  to  store  the 

program  and  the  system  information  are  illustrated  below. 


Each  magnetic  card  is  labeled  according  to  the  information 
stored  on  it.  In  the  upper  corners  of  each  card,  the  number 
in  each  space  is  to  indicate  the  bank  numbers  recorded  on  that 
card.  The  arrow  in  each  space  shows  which  direction  the  card 
must  be  inserted  into  the  calculator  when  reading  the  indicated 
bank.  The  space  across  the  center  of  the  card  is  available 
for  the  program  title  and  other  pertinent  information  such  as 
the  required  partition.  Below  this  line  are  two  rows  of 
boxes.  The  bottom  five  boxes  are  used  to  indicate  the  function 
of  the  defined  keys  £aJ  through  E  within  the  recorded  program. 
The  upper  row  of  boxes  is  used  similarly  for  the  keys  A' 
through  . 
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into 


To  read  a  magnetic  card,  press  CLRj  and  insert  the  card 


the  lower  slot  on  the  right  side  of  the  calculator.  The  drive 
motor  of  the  calculator  will  automatically  pull  the  card 
through  the  calculator.  The  number  of  the  bank  recorded  from 
the  entered  card  is  shown  in  the  display  after  that  card  side 
has  been  read  successfully. 


Repeat  the  procedure  (Press  | CLR  |  and  enter  a  card  side)  until 


all  four  banks  (1,  2,  3,  and  4)  have  been  correctly  read  (in 
any  order). 


Up  to  now,  the  computer  program  and  the  system  information  have 
been  read  into  the  calculator  from  the  magnetic  cards  and  will  be  stored 
in  the  calculator  as  long  as  the  switches  of  the  calculator  and  printer 
are  left  at  "ON”  positions  continuously. 


3  Enter  numerical  value  for  as  selected  generator  output  of 

TG3  and  press  .  Input  will  be  stored  into  data  register  20 
and  printed  out.  The  decrement  of  input  for  iteration  will  be 
in  the  display  after  a  very  short  period  of  computation. 
(During  computation,  a  "C"  appears  at  the  left  side  of  the 
otherwise  blank  display.) 


4 


Enter  the  value  of 
into  data  register 
will  appear  in  the 
tation. 


B 


Input  will  be  stored 


and  press 

21  and  printed  out.  The  decrement  of  input 
display  after  a  very  short  period  of  compu- 


5  Enter  the  value  of  W,.  and  press  [c] .  Input  will  be  stored 

into  data  register  22  and  printed  out.  The  decrement  of  input 
will  appear  in  the  display  after  a  very  short  period  of  compu¬ 
tation. 


Step 


6  Enter  the  value  of  Mu  and  press  D  .  Input  in  the  display  is 

n  — 

stored  into  data  register  23  and  printed  out. 


Enter  the  value  of  and  press  [_Ej .  Input  in  the  display  is 
stored  into  data  register  24  and  printed  out. 

Enter  the  value  of  D  and  press  2nd]  (~a|  .  Input  in  the  display 
is  stored  into  data  register  25  and  printed  out. 

Enter  the  value  of  and  press  2nd|  |~B~|  .  Input  the 
display  is  stored  into  data  register  26  and  printed  ~Sx^t . 

Enter  the  value  of  and  press  2nd j  ( C | .  Input  in  the 
display  is  stored  into  data  register  27  and  printed  out. 

Enter  the  value  of  P^.  and  press  2nd |  |  D | .  Input  in  the 
display  is  stored  into  data  register  28  and  printed  out. 

Enter  the  value  of  T  and  press  2nd|  |e| .  Input  in  the  display 
is  stored  into  data  register  29  and  printed  out. 

Note:  If  a  mistake  is  made  at  any  above  numerical  data  entering 

step  (3  to  12)  simply  press  CLR  and  start  over  again  for  that 
step . 
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To  run  the  program,  press: 

RST 

,  R/S 

When  the  program  is  in  execution,  there  is  a  blinking  "C"  on  the 
left  side  of  the  otherwise  blank  display.  The  following  results  will  be 
printed,  successively,  for  the  first  and  last  iterations: 


•  TG3  electrical  output,  kW 

•  TG4  electrical  output,  kW 

•  TG5  electrical  output,  kW 

•  Total  condensing  flow,  lb/hr 

•  Heat  rate,  Btu/kW-hr 

•  Unit  cost  of  self-generation,  C/kW-hr 

•  Total  cost  of  self-generation,  $ 

•  Unit  cost  of  purchased  power,  C/kW-hr 

•  Total  cost  of  purchased  power,  $ 

•  Total  cost,  $ 

•  Total  unit  cost,  C/kW-hr 


After  the  last  group  of  output  is  printed,  the  display  should  read 
"0"  indicating  that  the  program  has  been  successfully  run  and  has  stopped. 
A  sample  printout  is  given  in  Table  6.  The  first  outputs  are  the  results 
calculated  using  the  input  values  of  generator  outputs,  W^,  W^,  and 
given  by  the  user,  while  the  last  outputs  are  the  results  using  the 
optimum  self-generation  solutions  of  W^,  W^,  and  W^. 

For  repeated  or  subsequent  runs,  steps  3,  4,  and  5  (for  entering 
the  values  of  W^,  W^,  and  W,.)  are  required  to  be  repeated  whether  or  not 
the  entering  values  are  different  from  the  previous  values!  However, 
only  those  steps  from  6  to  12  need  to  be  repeated  for  which  entering 
values  differ  from  the  previous  run.  The  program  is  then  restarted  by 
pressing  RST  R/S  after  the  last  value  changed. 

Caution:  If  the  user  has  assumed  that  any  one  of  the  turbine/ 

generator  sets  has  been  shut  down  (by  entering  zeros  for  W^,  W^,  or  W,.) 
in  running  the  program,  he  must  rea*d  bank  no.  4  of  the  magnetic  card 
once  again  (as  described  in  step  2)  in  order  to  restore  the  system 
information  before  he  starts  to  rerun  the  program  for  a  different  combi¬ 
nation  of  turbine/generator  sets. 


CONCLUSIONS 


A  computer  program  designed  for  use  on  a  hand-held  TI-59  calculator, 
with  a  PC-100  printer,  has  been  developed  to  provide  the  personnel  of 
the  power  plant  at  SUBASE,  New  London,  Conn.,  with  a  rapid  computation 
of  the  optimum  operating  settings  for  the  turbine/generator  sets  during 
any  given  shift.  Input  parameters  of  this  computer  program  are  left 
free  and,  for  most  of  them,  are  not  required  to  be  entered  for  subsequent 
runs  unless  they  are  different  from  those  previously  entered.  The  most 
unique  feature  of  this  computer  program  is  that  it  provides  the  user 
with  maximum  flexibility  in  the  choice  of  the  mode  of  operation  of  the 
power  plant.  This  means  that  all  three  turbine/generator  sets  may  be 
assumed  online  for  operation,  or  any  other  combinations  of  one  or  two 
turbine/generators  may  be  assumed  online  in  the  process  of  determining 
the  optimum  mix  of  self-generated  and  purchased  electrical  power. 

The  economics  of  power  generation  at  SUBASE,  New  London,  Conn.,  has 
changed  quite  drastically  over  the  last  few  years  due  to  the  great 
increase  in  fuel  oil  costs.  Since  purchased  power  rates  have  not  increased 
in  the  same  proportion,  purchased  power  has  become  more  economical  than 
onsite  power  generation.  The  dramatic  influence  of  the  rate  structure 
and  fuel  price  on  the  total  unit  cost  is  evident  in  Figure  6.  The  data 
selected  for  use  in  Figure  6  are  shown  in  Table  7.  As  indicated  in 
Figure  6,  it  would  be  more  advantageous  to  the  SUBASE  to  purchase  elec¬ 
trical  power  from  the  utility  than  to  generate  it  onsite  while  the 
prices  of  fuel  and  utilities  are  increasing  at  the  same  rates  as  before. 
However,  due  to  the  very  nature  and  functioning  of  the  SUBASE,  the  power 
supply  reliability  dictates  a  reasonable  amount  of  base  generation  and 
spinning  reserve.  Therefore,  generation  assignments  should  be  made  in 
such  a  manner  as  to  strike  a  balance  between  reliability  requirements 
and  economics. 

The  significance  in  savings  resulting  from  the  optimum  operation  of 
the  existing  cogeneration  system  at  the  SUBASE  has  also  been  changed 
over  the  years.  Figure  6  shows  that  the  operation  of  the  power  plant  at 
the  optimum  make/buy  ratio  in  1972  was  not  as  critical  as  it  became  in 
1976  and  1979.  For  example,  curve  no.  1  illustrates  that  the  optimum 
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operating  point  in  1972  was  at  a  make/buy  ratio  3,500  kW/1 1,500  kW  for  a 
given  demand,  15,000  kW.  Even  when  operating  100%  off  the  optimum 
point,  i.e.,  7,000  kW/8,000  kW,  in  1972  there  would  have  resulted  an 
additional  total  unit  cost  of  about  0.06  cents  per  kW-hr.  However,  a 
similar  deviation  from  the  optimum  points  (100%  off)  for  curves  2,  3  and 
4  would  result  in  additional  total  unit  cost  of  0.30,  0.40  and  0.65 
cents  per  kW-hr  respectively,  the  latter  of  which  is  typical  of  1979. 

Hence,  it  can  be  concluded  that  more  savings  would  be  achieved  by  operating 
the  existing  power  plant  at,  or  close  to,  the  determined  optimum  make/buy 
decision  points  as  the  fuel  price  is  expected  to  increase  continuously. 

Figures  7  and  8  present  the  results  of  the  parametric  study  of  the 
effect  of  steam  loads  and  fuel  price  on  total  unit  cost,  respectively. 

As  clearly  seen,  the  variations  of  steam  demand  will  affect  the  location 
of  minimum  points  and  the  total  unit  cost  levels  dramatically.  However, 
the  increase  of  fuel  prices  may  not  affect  the  locations  of  minimum 
points  but  would  result  in  higher  additional  total  unit  costs  if  the 
plant  was  not  operated  at  the  minimum  point. 

RECOMMENDATIONS 

In  order  to  evaluate  the  savings  resulting  from  the  implementation 
of  this  program,  a  user  data  sheet  is  designed  and  shown  in  Table  8  for 
users  to  record  all  the  relevant  information.  There  are  four  columns  in 
this  user  data  sheet.  The  first  column  records  the  date  and  time  when 
the  operator  starts  to  use  this  program.  The  second  column  is  used  to 
record  the  operational  information.  The  third  column  is  for  the  results 
printed  out  from  the  printer.  The  last  column,  which  is  very  important 
for  the  evaluation,  is  for  the  actual  information,  after  the  fact,  that 
will  be  used  to  calculate  the  costs  to  be  compared  with  those  in  the 
third  column. 

Although  the  computer  program  developed  in  this  study  is  for  use  at 
SUBASE,  New  London,  Conn.,  its  applicability  to  other  Navy  activities 
which  already  practice  cogeneration  with  steam  turbine/generator  equipment 


is  recommended  to  be  investigated  by  modifying  this  program  according  to 
the  site-specific  requirements,  such  as  thermodynamic  characteristics  of 
the  cogeneration  system,  the  utility  rate  structure  and  load  demands. 


Table  2.  Utility  Structure  at  SUBASE,  New  Longon,  Conn, 
(as  of  April  1980) 


Demand  Charge  per  Month: 

(with  . 15  per  kW  credit) 

$4. 72/kW 

Energy  Charge: 

(subject  to  5%  discount) 

First  50,000  kW-hr 

$0 . 0343/kW-hr 

Over  50,000  kW-hr  and  up 
to  300  x  billing  demand 

$0 . 0269/kW-hr 

Over  300  x  billing  demand 

$0.021 1/kW-hr 
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Table  i.  ConLinued 
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Table  4.  System  Information  (Refer  to  Nomenclature 
List  for  Designations) 


Designations 

Value 

Data 

Registers 

Units 

Hv 

150,000 

00 

Btu/gal 

hs”he 

435 

01 

Btu/lb 

hs:hh 

100 

02 

Btu/lb 

Vhl 

205 

03 

Btu/lb 

he‘hc 

860 

04 

Btu/lb 

E3T 

0.70 

05 

X 

E4T 

0.68 

06 

X 

E5T 

0.72 

07 

% 

M3M 

2,000 

08 

lb/hr 

M4M 

2,300 

09 

lb/hr 

M5M 

5,000 

10 

lb/hr 

W3M 

1,500 

11 

kW 

W4M 

2,000 

12 

kW 

W5M 

1,500 

13 

kW 

eb 

0.72 

14 

% 

3,413 

15 

conversion  factor 

Du 

12,654 

16 

kW 

0.01 

17 

factor 

Table  5.  User  Instructions 


Parameters 

W3  (A) 

W4  (B) 
w5  (C) 

Mh  (D) 

ml  (e) 

D  (A') 

Kj  (B’) 

fa  (O 
PF  (D ' ) 

T  (E’) 


Input 

3,500 

4,000 

5,000 

81,300 

21,600 

16,500 

100,000 

0.0024 

58 

8 

First  Output 

3,500 

4,000 

5,000 

110,999.68 

12,504.21 

6.72 

6,715.22 

5.01 

1,601.68 

8,316.90 

6.30 

Final  Output 

2,100 

2,400 

3,000 

55,052.22 

11,180.08 


continued 


Table  6.  Continued 


Parameters 


C 

C 

C 

C 

C 

C 


US 

TS 

UP 

TP 

T 

U 
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Final  Output 

6.00 

3,602.47 

3.85 

2,773.68 

6,376.15 

4.83 


Table  7.  Utility  Rate  Structures  and  Fuel  Prices  Used  in  Figure  6 


Curve 

#1 

Curve 

#2 

Curve 

#3 

Curve 

#4 

1972 

March 

1976 

Aug  1979 

Before 

After 

Demand  charge,  $/kW 
(with  discount) 

B 

2.475 

4.15 

4.72 

Energy  charge,  C/kW-hr 

(subject  to  5%  discount) 

B 

1st  50,000 

Ifls 

1.95 

3.19 

3.43 

over  50,000  and 
up  to  300  x  demand 

1.02 

1.53 

2.50 

2.69 

over  300  x  demand 

0.8 

1.2 

1.96 

2.11 

fuel  price,  C/gal 

10 

32 

32 

58 

34 


item  ;it  SUBASK,  New  London,  Conn, 
symbols) 


maximum  throttle  flow 


igure  2.  Turbine/generator  No.  3  (TC.3)  performance. 


Throttle  Flow  <103  Ib/hr) 


Totsl  Unit  Cost  (eents/kW-hr) 


Total  Unit  Co«t  (centt/kW-hr) 


curve  1 


Self-Generation  to  Meet  Demand  of  15,000  kW  (kW) 


NOMENCLATURE  LIST 


c 

Total  cost  ($) 

HV 

Heating  value  of  fuel 

1 

V 

(Btu/gal ) 

CTP 

Total  cost  ($)  of  purchased 
power 

Kp 

Purchased  electrical  power 

r 

(kW) 

CTS 

Total  cost  ($)  of  self¬ 
generated  power 

Ki 

Cumulative  purchased  energy 
(kW-hr)  at  the  beginning  of 

cu 

Unit  cost  (C/kW-hr) 

each  run 

CUP 

Unit  cost  (C/kW-hr)  of 

K2 

Cumulative  purchased  energy 

purchased  power 

(kW-hr)  at  the  end  of  each 

run 

cus 

Unit  cost  (C/kW-hr)  of  self¬ 
generated  power 

Mc 

Total  condensing  flow  from 

turbine  (lb/hr) 

D 

Electrical  demand  (kW) 

mh 

High-pressure  (200-psig) 

Du 

Utility  peak  demand  (kW) 

export  steam  (lb/hr) 

eb 

Boiler  overall  efficiency 

MiE 

Flow  rate  (lb/hr)  at  turbine 

outlet  of  TGi 

EiT 

Overall  efficiency  of  TGi 

MiH 

Flow  rate  (lb/hr)  at  high- 

fa 

Fuel  adjustment  factor 

pressure  (200-psig) 

($/kW-hr ) 

extraction  point  of  TGi 

Hc 

Enthalpy  (Btu/lb)  of  the 

Mil 

Flow  rate  (lb/hr)  at  low- 

mixed  condensate 

pressure  (5-psig) 
extraction  point  of  TGi 

Enthalpy  (Btu/lb)  of  the 
mixed  exhaust  steam 

MiM 

Minimum  required  exhaust 
steam  flow  (lb/hr)  of  TGi 

HiC 

Enthalpy  (Btu/lb)  at 
condenser  outlet  of  TGi 

MiT 

Flow  rate  (lb/hr)  at  turbine 
throttle  of  TGi 

HiE 

Enthalpy  (Btu/lb)  at  turbine 
outlet  of  TGi 

ml 

Low-pressure  (5-psig)  export 

steam  (lb/hr) 

«H 

Enthalpy  (Btu/lb)  of  high- 
pressure  (200-psig)  export 

PF 

Fuel  price  (C/gal) 

HiH 

steam 

Enthalpy  (Btu/lb)  at  high- 
pressure  (200-psig) 
extraction  point  of  TGi 

T 

Time  period  of  interest  (hr) 
for  each  run  which  may  be  a 
period  of  one  shift  or  less 

T 

HR 

Heat  rate  of  turbine/ 

hl 

Enthalpy  (Btu/lb)  of  low- 

generator  sets 

la 

pressure  (5-psig)  export 
steam 

W. 

l 

Electrical  output  (kW)  of  TG: 

H.t 

iL 

Enthalpy  (Btu/lb)  at  low- 
pressure  (5-psig)  extraction 
point  of  TGi 

W.„ 

iM 

Low  limit  of  electrical 
output  (kW)  of  TGi 

WT 

Total  self-generated 

HS 

Enthalpy  (Btu/lb)  of  steam 

electricity  (kW) 

from  boilers 

H Enthalpy  (Btu/lb)  at  turbine 
11  throttle  of  TGi 

^/Y6 
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Great  l  akes  It:  l  ode  10.  Great  Lakes.  II  Code  110.  Oakland.  CA:  Code  120.  Oakland  CA:  Code  1201  . 

( Librarv)  San  Diego.  CA;  Code  154.  Great  Lakes.  II  :  Code  2tHi.  Great  Lakes  It:  Code  220  Oakland.  CA; 
Code  220.1.  Norlolk  VA:  l 'ode  30C.  San  Diego.  CA;  I'ode  4IHI.  Great  Lakes.  II.;  Code  4(H).  Oakland.  CA: 
Code  400.  Pearl  Harbor.  HI;  Code  400.  San  Diego.  CA:  Code  420.  Great  Lakes.  IL;  l  ode  420.  Oakland. 
CA;  lode  42B  (R.  Paseua).  Pearl  Harbor  HI;  Code  505A  (11  Wheeler);  Code  600,  Great  Lakes.  IL:  Code 
601.  Oakland.  CA:  Code  610.  San  Diego  l  a.  Code  7(10.  Great  Lakes.  IL:  I.FJG  I  I  MeClaine.  Yokosuka. 
Japan:  Utilities  Officer.  Guam:  NO  (Code  201  Oakland.  CA 
TV  A  Smelser.  Knoxville,  l'enn. 

C  S  MERCHANT  MARINE  ACADEMY  Kings  Point.  NY  (Reprint  Custodian) 

Lis  GEOLOGICAL  SURVEY  Oft  Marine  Geology.  Piteleki.  Reston  VA 

USAF  Jack  S.  Speneer.  Washington.  DC 

USAF  REGIONAL  HOSPITAL  Fairchild  AFB.  WA 

CSCG  (G-ECV)  Washington  D e:  (Smith).  Washington.  DC 

IISCG  RAD  CENTER  Tech.  Dir.  Groton.'  IT 

USNA  Ch  Meeh  Engr  Dept  Annapolis  MD;  Energy-Eniiron  Study  Grp.  Annapolis.  Ml);  tngr.  Dii.  (C.  Wu) 
Annaplolis  Ml):  Environ  Prot.  RAD  Prog  (J.  Williams).  Annapolis  MD.  Ocean  Sys.  Eng  Dept  (Dr 
Monney)  Annapolis.  MD.  PWD  Engr.  Div.  (C.  Bradford)  Annapolis  MD 
ARIZONA  State  Energy  Programs  Off  .  Phoenix  AZ 

BONNEVILLE  POWER  ADMIN  Portland  OR  (Energy  Consrv  Off..  D.  Davev) 

BROOKIIAVEN  NAIL  LAB  M  Steinberg.  Upton  NY 

CALIFORNIA  STATE  UNIVERSITY  LONG  BEACH.  I  A  (CHEI.APA  1 1) 

CORNELL  UNIVERSITY  Ithaca  NY  (Serials  Dept.  Engr  Lib.) 

DAMES  A  MOORE  LIBRARY  LOS  ANGELES.  CA 
FLORIDA  ATLANTIC  UNIVERSITY  Boca  Raton.  FI  (McAllister) 

FLORIDA  TECHNOLOGICAL  UNIVERSITY  ORLANDO.  FI.  (HARTMAN) 

FOREST  INST.  FOR  OCEAN  A  MOUNTAIN  Carson  City  NV  (Studies  -  Library) 

FUEL  A  ENERGY  OFFICE  CHARLESTON.  WV 

GEORGIA  INSTITUTE  OF  TECHNOLOGY  (IT  R  Johnson)  Atlanta.  GA 
HAWAII  STATE  DEPI  OF  PLAN.  A  ECON  DEV.  Honolulu  HI  (Tech  Info  Ctr) 

INDIANA  ENERGY  OFFICE  Energy  Group.  Indianapolis.  IN 
WOODS  HOLE  OCEANOGRAPHIC  INST  Woods  Hole  MA  (Winget) 


M  I  NI  SIAM  COLLEGE  Kccni:  Nil  (C  unnmgh.iml 

I  l  mull  UNIVERSITY  HI- 1  III  l  III  SI.  PA  (MARIM  (.lOlHIINU  \l  I  \B  RU  II  \KI  )S ) .  Bethlehem 
PA  (l.tndcrmun  l  ib  No  .'ll.  Ilciksteiner  I 
MAINE  OFFICE  OK  I  NI  RGY  RESOURCES  Augusta.  Ml 
MISSOURI  ENERGY  AGENCY  Icltcrson  Cm  MO 

Mil  Cambridge  MA  iRm  III*  SOU.  lech  Reports.  I  net  lib  i  I  .imbndge  MA  iHurlem.ini 

MONTANA  ENERGY  OFFICE  Anderson.  Helena.  Ml 

NA  I  I  ACADEMY  OK  l  NG  Al  I  XANDRIA.  \  A  (SI  \KI  I  IK  i 

NEW  HAMPSHIRE  Concord.  Nil.  (Uoicrnor  s  Council  On  Inctgyi 

NEW  MEXICO  SOLAR  I  NI  ROY  INS  I  l)r  /wibcl  I  .is  (  rims  N\| 

NY  CITY  COMMUNITY  COI  .LEGE  BROOKLYN.  NY  (IIBRARSi 
NYS  ENERGY  OEFICE  Library.  Albany  N\ 

POLLUTION  ABATEMENT  ASSOC  Graham 
PURDUE  UNIVERSITY  Lafayette.  IN  ICE  Engr  l  ib) 

CONNECTICUT  Hartford  Cl  (Dept  of  Plan  &  I  nergv  Police) 

SCRIPPS  INSTITUTE  OF  OCEANOGRAPHY  LA  JOLLA.  CA  (ADAMSl 

SEATTLE  l'  Prof  Schwaegler  Seattle  WA 

STANFORD  UNIVERSITY  Engr  Lib.  Stanlord  CA 

STATE  I'NIV  OF  NEW  YORK  Fort  Schuyler.  NY  (tongohardil 

LEX  AS  Ait  M  UNIVERSITY  W.B.  Ledbetter  College  Station.  IX 

UNIVERSITY  OF  C  ALIFORNIA  Energy  Engineer.  Davis  CA;  LIVERMORE.  CA  | LAWRENCE 
LIVERMORE  LAB.  TOKAR/.I 

UNIVERSITY  OF  DELAWARE  Newark.  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  HONOLULU.  HI  (SCIENCE  AND  LEC  H  DIV  I 
UNIVERSITY  OF  ILLINOIS  URBAN  A.  IL  (LIBRARY):  URBAN  A.  II  (NEWMARK) 

UNIVERSITY  OF  MASSACHUSETTS  (Heronemus).  Amherst  MA  C  E  Dept 
UNIVERSITY  OF  NEBRASKA-LINCOLN  Lincoln.  NE  (Ross  Ice  Shelf  Prop  I 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Set  (Library).  Port  Arkansas  LX 
UNIVERSITY  OF  TEXAS  AT  AUSTIN  AUSTIN.  LX  (THOMPSON) 

UNIVERSITY  OF  WASHINGTON  Seattle  WA  (E.  Linger);  Seattle.  WA  Transportation.  Construction  X 
Geom.  Div 

IfNIVERSITY  OF  WISCONSIN  Milwaukee  WI  (Ctr  of  CJrcat  Lakes  Studies) 

VIRGINIA  INST.  OF  MARINE  SCI.  Gloucester  Point  VA  (Library) 

ARVID  GRANT  OLYMPIA.  WA 

ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMITH) 

BAGGS  ASSOC.  Beaufort.  SC 

BECHTEL  CORP.  SAN  FRA  NCI  SCO.  CA  (PHELPS) 

BRITISH  EMBASSY  Sei.  &  Tech.  Dept.  (J.  McAuley).  Washington  DC 
BROWN  &  CALDWELL.  E  M  Saunders  Walnut  Creek.  CA 
BROWN  &  ROOT  Houston  TX  (D.  Ward) 

CANADA  Nova  Scotia  Rsch  Found.  Corp  Dartmouth.  Nova  Scotia:  Trans  Mnt  Oil  Pipe  Lone  Corp. 
Vancouver.  BC  Canada 

CTTEMED  CORP  Lake  Zurich  IL  (Dearborn  Chcm  Div.I.ib.) 

COLUMBIA  GULF  TRANSMISSION  CO  HOUSTON.  TX  (ENCi.  LIB  ) 

DESIGN  SERVICES  Beck.  Ventura.  CA 
DIXIE  DIVING  CENTER  Decatur.  CIA 
DURLACH.  O'NEAL.  JENKINS  &  ASSOC  Columbia  SC 
FORD.  BACON  &  DAVIS.  INC.  New  York  (Library) 

MCDONNEL  AIRCRAFT  CO.  Dept  501  (RTF  Fayman).  St  Louts  MO 
MEDERMOTT  &  CO.  Diving  Division.  Harvey.  LA 

NEWPORT  NEWS  SHIPBLDG  &  DRYDOCK  CO  Newport  News  VA  (  lech  Lib  ) 

PACIFIC  MARINE  TECHNOLOGY  Duvall.  WA  (Wagner) 

PORTLAND  CEMENT  ASSOC,  Skokie  IL  (Rsch  &  Dev  Lab.  Lib.) 

RAYMOND  INTERNATIONAL  INC  E  Colic  Soil  Tech  Dept.  Pennsauken.  NJ 
SANDIA  LABORATORIES  Albuquerque.  NM  (Vortman):  Library  Div  .  Livermore  CA 
SCHUPACK  ASSOC  SO  NORWALK.  CT  (SCHUPAC  K) 

SEATECH  CORP  MIAMI.  FL  (PERONI) 

SHELL  DEVELOPMENT  CO.  Houston  TX  (C  Sellars  Jr  ) 
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SWI  1)1  N  VBB  (librarx).  Stockholm 

I  EX  IRON  INC  Bl  FEALO.  MV  (R!  SI  ARC  II  (T  Nil  R  [.IB.) 

I  RW  S')  SI  I  MS  REDONDO  BEACH.  CA  (DAI) 

CM  I  ED  KINODOM  Cement  A  Concrete  Assoc  Wcxh.im  Springs.  Slough  Bucks;  Library.  Bristol:  R  Browne. 
Southall.  Middlesex .  lax  lor.  Woodrow  C'onstr  (0I4PI.  Southall.  Middlesex.  Taylor.  Woodrow  C'onxtr 
(Smith).  Southall.  Middlesex;  l  nix  of  Bristol  (R  Morgan).  Bristol 
CM  I  ED  IT  CHNOI  OC.IES  W  indsor  Locks  C  I  (Hamilton  Std  Dix  .  I.ihrarv) 

W  ARD.  WOISTTNHOI  D  ARC  HITECTS  Sacramento.  C  A 

WESITNCiHOl  SE  ELECTRIC  CORE  Annapolis  MD  (Oceanic  Dix  lab.  Brx.inl:  l.ihrarx.  Pittshurgh  PA 
WM  CLAPP  LABS  HAITI  I  I  I  DCXBIRY.  MA  (LIBRARY) 

BRAHTZ  La  Jolla.  CA 
KE  IRON.  BOB  Ft  Worth.  TX 
KRl  /.IC  .  I  P  Silver  Spring.  MD 
C  APT  MCRPIIY  Sunnyvale.  CA 
T  W  MERMEL  Washington  DC 
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